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Abstract 

We show that all nontrivial CP-violating phases can be determined in terms 
of three lepton flavor mixing angles and the ratio of Am^ un to Am^ tm in 
the minimal seesaw model in which the Frampton-Glashow-Yanagida (FGY) 
ansatz is incorporated. This important point allows us to make very specific 
predictions for the cosmological baryon number asymmetry and CP violation 
in neutrino oscillations. A measurement of the smallest neutrino mixing angle 
will sensitively test the FGY ansatz, in particular in the case that three light 
neutrinos have a normal mass hierarchy. 
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I. INTRODUCTION 



In the minimal standard model of electroweak interactions, the lepton number conser- 
vation is assumed and neutrinos are exactly massless Weyl particles. However, today's 
Super-Kamiokande [1], SNO [2], KamLAND [3] and K2K [4] neutrino oscillation experi- 
ments have provided us with very strong evidence that neutrinos are actually massive and 
lepton flavor mixing does exist. The most economical modification of the minimal standard 
model, which can both accommodate neutrino masses and allow lepton number violation to 
explain the cosmological baryon asymmetry via leptogenesis [5], is to introduce two heavy 
right-handed neutrinos N 1>2 and keep the Lagrangian of electroweak interactions invariant 
under the SU(2)l x U(1)y gauge transformation [6,7]. After the spontaneous electroweak 
symmetry breaking, this simple but interesting model yields the following neutrino mass 
term: 

-£ mass = {v^v T ) M D + \{Nl,N$) Mr (^) + h.c. , (1) 

where N£ = CNf with C being the charge-conjugation operator; and (u e , z/ M , u T ) denote the 
left-handed neutrinos. The Dirac neutrino mass matrix is a 3 x 2 rectangular matrix, 
and the Majorana neutrino mass matrix Mr is a 2 x 2 symmetric matrix. The scale of 
M D is characterized by the electroweak scale v — 174 GeV. In contrast, the scale of M R 
can be much higher than v, because Ni and N 2 are SU(2) L singlets and their corresponding 
mass term is not subject to the scale of gauge symmetry breaking. Then one may obtain 
the effective (light and left-handed) neutrino mass matrix M v via the well-known seesaw 
mechanism [8]: 

M v « M d MrX • ( 2 ) 

Without loss of generality, both Mr and the charged lepton mass matrix M\ can be taken 
to be diagonal, real and positive; i.e., M R = Diag{M 1 ,M 2 } and M = Diag{m e , m M , m T }, 
where M 12 denote the masses of two heavy Majorana neutrinos. In such a specific flavor 
basis, the low-energy neutrino phenomenology is governed by M u , while the cosmological 
baryon number asymmetry is associated with Mq via the leptogenesis mechanism. 

Unfortunately, the model itself has no restriction on the structure of Md. In Ref. [6], 
Frampton, Glashow and Yanagida (FGY) have conjectured that M D may take the form 
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The texture zeros in Md could stem from an underlying horizontal flavor symmetry. With 
the help of Eq. (2), one may straightforwardly arrive at 
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from Eq. (3); or 
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from Eq. (4). Note that Det(M„) = holds in either case 1 . Note also that \Det(M u )\ = 
m 1 m 2 m3 holds in the chosen flavor basis, where m 8 (for % = 1,2, 3) denote the masses of 
three light neutrinos. Thus one of three neutrino masses must vanish. As the solar neutrino 
oscillation data have set m 2 > mi [2], we are left with either mi = (normal hierarchy) 
or m 3 = (inverted hierarchy). In Ref. [9], the Majorana neutrino mass matrix with one 
texture zero and one vanishing eigenvalue has been classified and discussed in some detail. 

The main purpose of this paper is to reveal a very striking feature of the minimal seesaw 
model in which the FGY ansatz is incorporated: all nontrivial CP-violating phases can be 
calculated in terms of the lepton flavor mixing angles (9 x ,9 y ,9 z ) and the ratio of Am 2 un to 
Am^ tm , where Am^ and Am^ tm stand respectively for the typical mass-sqaured differences 
of solar and atmospheric neutrino oscillations. This important point, which was not observed 
in the previous analyses of the minimal seesaw model [6,7,10], implies that a stringent test 
of the FGY ansatz can simply be realized once the smallest mixing angle 6 Z is measured 
or constrained to a better degree of accuracy. Considering both normal and inverted mass 
hierarchies of three light neutrinos, we obtain very specific predictions for the cosmological 
baryon number asymmetry, the effective mass of neutrinoless double beta decay and CP 
violation in neutrino oscillations. 



II. DETERMINATION OF CP- VIOLATING PHASES 

Without loss of generality, one may always redefine the phases of charged lepton fields 
to make a, b and b' of M D real and positive [6]. In other words, only a' is complex and its 
phase 4> = arg(a') is the sole source of CP violation in the model under discussion. Because 



x It has been shown in Ref. [9] that Det(M^) = is independent of the specific texture zeros taken 
in Eq. (3) or (4). In other words, the minimal seesaw model itself guarantees that Det(M^) = 
holds automatically. 
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both Mi and M R have been taken to be diagonal, real and positive, M v can in general be 
parametrized as follows: 



M v = (PUQ) 



fm 1 \ 

m 2 
V mj 



(PUQf , (7) 



where P = Diag{e iCt , e t/3 , e n } and Q = Diag{e* p , e ta , e tu } are two phase matrices, and U is 
given by 



U 



f c x c z s x c z s z 

+ C X Cy6 SyC Z 

(8) 

\ CxCy^z ~t~ Sx^y^ Sx^y^z Cx^y^ CyCz 



with s x = sin^, c x = cos^ and so on. The phase parameter of U (Dirac phase) governs the 
strength of CP violation in neutrino oscillations, while two independent phase parameters 
of Q (Majorana phases) are relevant to the neutrinoless double beta decay [11]. The phases 
of P cannot be neglected in the parametrization of M v — their essential role is to fulfil a 
complete match between the phases of M v in Eq. (5) or (6) and those defined in Eqs. (7) 
and (8). It is then obvious that the nontrivial phases of P, U and Q should have definite 
relations with 0. 

Note that three mixing angles of U can directly be given in terms of the mixing angles 
of solar, atmospheric and reactor [12] neutrino oscillations. Namely, 9 X pa 6 sun , 9 y pa # atm 
and 9 Z ~ # c h z hold as a good approximation. In view of current experimental data, we have 
9 X pa 32° and y pa 45° (best-fit values [13]) as well as 9 Z < 12°. The mass-squared differences 
of solar and atmospheric neutrino oscillations are defined respectively as Am^ = m 2 , — m\ 
and Am^ tm = \m\ — m||. Their best-fit values read Am 2 un pa 7.13 x 10 -5 eV 2 and Am 2 tm pa 
2.6 x 10~ 3 eV 2 [13]. These typical numbers will be used in our numerical calculations. 

To be more specific, we shall concentrate on the FGY ansatz for Md in Eq. (3) or 
equivalently M v in Eq. (5). Some brief comments will be given on the consequences of 
in Eq. (4) or equivalently M v in Eq. (6). Indeed, both possibilities lead to very similar 
phenomenological results. 



A. Normal neutrino mass hierarchy (m\ = 0) 



If mi = holds, we obtain m 2 = y Am 2 un pa 8.4 x 10 3 eV and m 3 = y Am 2 un + Am 2 tm ~ 
5.2 x 10~ 2 eV. In this case, only a single Majorana phase of CP violation is physically 
nontrivial. Hence the phase matrix Q can be simplied to Q — Diag{l, e lcr , 1} and six 
independent matrix elements of M v can be written as 



(M„) n 


= e 2ia 


2 2 2ia i 2~ 

m 2 s x c z e + m 3 s z 


j 






(M„) 22 


= e 2l(i 


Wl2 ^ Sx^y^z C x C y € 




(M„) 33 


= e 2n 


TY12 {^xCySz ~t~ C x S y C 




2 e 2la + m 3 c 2 c 2 





4 



and 

(M„) 12 = e^ a+ V 
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Because of (M u )i 3 = as shown in Eq. (5), we straightforwardly obtain 
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where £ = m 2 /m 3 0.16. This result implies that both 5 and ex can definitely be determined, 
if and only if the smallest mixing angle 9 Z is measured. 

To establish the relationship between <b and S, we need to figure out a, (3 and 7. Because 
a, 6 and 6' are real and positive, (M u )u, (M u ) 23 and (M v ) 33 must be real and positive. Then 
a, /3 and 7 can be derived from Eqs. (10) and (11): 
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Then the overall phase of (M u ) 12 , which is equal to the phase of a', is given by 



= a + (3 — arctan 
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Now that all six phase parameters (5, a, <b, a, (3 and 7) have been determined in terms of 
£, 9 X , 9 y and 9 Z , a measurement of the unknown angle 9 Z becomes crucial to test the model. 



B. Inverted neutrino mass hierarchy (m 3 = 0) 

If m 3 — holds, we arrive at mi = \J Am^ — Am 2 un ps 5.0 x 10~ 2 eV and m 2 



atm 



5.1 x 10~ 2 eV. Since only a single Majorana phase of CP violation is physically 
nontrivial, we can always simplify the phase matrix Q to Q = Diag{l, e" 7 , 1}. In this case, 
six independent matrix elements of M v turn out to be 
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and 
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where £ = m\jm2 ~ 0.98. Once the smallest mixing angle # z is observed, one may determine 
both 5 and a with the help of Eq. (16). 

The phase parameters a, (3 and 7 can be fixed by taking account of the positiveness of 
(M v )n, (M„) 23 and (M u ) 33 . With the help of Eqs. (14) and (15), we obtain 
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Then the overall phase of (M u )i 2 , which is equal to the phase of a', is given by 



= a + (3 — arctan 
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Again, a measurement of the unknown mixing angle 9 Z will allow us to determine all six 
phase parameters (<5, a, <f>, a, (3 and 7). 



C. Numerical dependence of (5, a, 0, a, (3, 7) on 2 

Using the best-fit values of Am^, Am^ tm , ^ and 6 y , we illustrate the numerical depen- 
dence of six phase parameters (<5, a, 0, a, (3, 7) on the smallest mixing angle # z in Fig. 1(a) 
and Fig. 1(b) for the mi = case and in Fig. 2(a) and Fig. 2(b) for the m 3 = case. Some 
discussions are in order. 

(1) In the mi = case, 9 Z is restricted to a very narrow range 4.0° < 9 Z < 4.4° (namely, 
0.070 < s z < 0.077). This result implies that the FGY ansatz with mi = is highly sensitive 
to 9 Z and can easily be ruled out if the experimental value of 9 Z does not really lie in the 
predicted region. To a good degree of accuracy, we obtain 5 2a, <fi ps a — <r, /3 p=j —7, 
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and 7^0. These instructive relations can essentially be observed from Eqs. (11), (12) and 
(13), because of s z <C 1. Note that we have only shown the dependence of 5 on Z in the 
range < 5 < ir. The reason is simply that only this range can lead to Y B > (i.e., the 
positive cosmological baryon number asymmetry), as one can see later on. 

(2) In the m 3 = case, there is no strong constraint on 9 Z except that 9 Z > 0.36° (or 
equivalently s z > 0.0063) must hold. We see that 5~/9~0 + 7r, 0^7, cr ~ —a and 
a k hold to a good degree of accuracy. These results can also be observed from Eqs. (16), 
(17) and (18) by taking account of s z <C 1. Again, we have used the positive sign of 1b 1° 
constrain the allowed range of 5. The other phase parameters are required to take possible 
values between — n and +7r. 

As useful by-products, the Jarlskog parameter of CP violation (Jcp [14]) and the effective 
mass of neutrinoless double beta decay ({m) ee [15]) can be calculated. We show the numerical 
results of Jcp versus (m) ee in Fig. 1(c) for mi = and in Fig. 2(c) for m 3 = 0, respectively. 
It is clear that < J C p ^ 0.016 and 2.1 meV < (m) ee < 2.7 meV in the m x = case, 

while < J C p ^ 0.035 and (m) ee « \J Aml tm « 0.05 eV in the m 3 = case. The present 
experimental upper bound of (m) ee is (m) ee < 0.35 eV at the 90% confidence level [16]. 



III. BARYON ASYMMETRY VIA LEPTOGENESIS 

Because of lepton number violation, two heavy Majorana neutrinos N t (for i = 1 and 2) 
may decay into IH and its CP-conjugate state, where / denotes the left-handed lepton doublet 
and H stands for the Higgs-boson weak isodoublet. The decay occurs at both the tree level 
and the one-loop level (via self-energy and vertex corrections), and their interference leads 
to a CP- violating asymmetry E{ between the CP-conjugated iVj — > I + H and Ni — > I + H* 
processes [5]. If the masses of Aq and N 2 are hierarchical (i.e., Mi M 2 ), the interactions 
of Ni can be in thermal equilibrium when decays. The asymmetry e 2 is therefore erased 
before Aq decays, and only the asymmetry e\ produced by the out-of-equilibrium decay of 
Ni survives. In the flavor basis chosen above, we have 



__ T(Ni -+1 + H) - T(Ni -+1 + H*) 
= T(Ni ^l + H) + Y(Ni 1 + H*) 

3 Mi Im[(M^M D ) 12 
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In deriving this formula, we have used (MdMd)h = a 2 + \a' | 2 and (M^M D )i 2 = (a')*b as 
well as 
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One can see that E\ is independent of M 2 , as long as M 2 ^> M\ is satisfied. A nonvanishing 
E\ may result in a net lepton number asymmetry Yl = n h/ s = d£i/g*, where = 106.75 
is an effective number characterizing the relativistic degrees of freedom which contribute to 
the entropy s of the early universe, and d accounts for the dilution effects induced by the 
lepton- number- violating wash-out processes [17]. If the effective neutrino mass parameter 
fhi = (Mjt,M D )ii/Mi [18] lies in the range 10~ 2 eV < m x < 10 3 eV, one may estimate the 
value of d by use of the following approximate formula [19] 2 : 

/l(T 3 eV\ r / mi \ 

d « 0.3 In -^—) 

\ mi J [ U0- 3 eV/ 

The lepton number asymmetry Y~l is eventually converted into a net baryon number asym- 
metry Y B via the nonperturbative sphaleron processes [21]: Y B = n B /s — 0.55Y l . A 
generous range 0.7 x 10~ 10 < F B ~ 1-0 x 10 -10 has been drawn from the recent WMAP 
observational data [22]. 

It is clear that E\ and Y B only involve two free parameters: Mi and <fi. Because is 
associated with the unknown flavor mixing angle 9 Z , one may analyze the dependence of Y B 
on 9 Z for given values of M\. For mi = and m 3 = cases, we plot the numerical results of 
1b in Fig. 1(d) and Fig. 2(d) respectively. Some comments are in order. 

(1) In the mi = case, current observational data of Y B require Mi > 2.9 x 10 10 GeV 
for the allowed ranges of s z . Once s z is precisely measured, it is possible to fix the value of 
Mi in most cases (e.g., Mi = 10 11 GeV will be ruled out, if s z ss 0.074 holds). 

(2) In the m 3 = case, Mi > 2.7 x 10 13 GeV is required by current observational data 
of Y B . Although 6 Z is less restricted in this scenario, it remains possible to determine the 
value of Mi once 9 Z is measured (e.g., Mi ps 5 x 10 13 GeV is expected, if s z m 0.1 holds). 

(3) One can carry out a similar analysis of Y B in the framework of supersymmetric seesaw 
and leptogenesis models. However, the FGY ansatz does not favor Mi < 10 8 GeV, which 
crucially affects the maximum reheating temperature of the universe after inflation in the 
generic supergravity models [17]. 

We remark that the cosmological baryon number asymmetry is closely correlated with 
the Jarlskog parameter of CP violation in the FGY ansatz 3 . For illustration, we plot the 
numerical correlation between Y B and Jcp in Fig. 3, where Mi = 5 x 10 10 GeV for the 
mi = case and Mi = 5 x 10 13 GeV for the m 3 = case have typically been taken. One 
can see that the observationally-allowed range of Y B corresponds to J C p ~ 1% in the m x = 



(21) 



2 For Mi 10 14 GeV, Giudice et al have presented a different approximate formula for the dilution 

r in 1 3.3 x 10" 3 eV / fhi \ 116 TT7 n , , 

factor (denoted as ri 20 ): - f» h -. — — . We find that the ratio n a 

r] mi \5.5xlO" 4 eV/ 

will vary in the range 0.9 < rj/d < 2.0, if rh\ takes values in the region 10 -2 eV < m\ < 10 3 eV. 
Thus there is no significant inconsistency between two empirical formulas. 

3 A similar point has been discussed by Endoh et al [7] in a more general way for the minimal seesaw 
model with leptogenesis. See, also, Ref. [23] for a possible seesaw bridge between leptogenesis and 
CP violation at low energies beyond the minimal seesaw model. 
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case and Jcp ~ 2% in the m 3 = case. The correlation between Y& and Jcp is so strong 
that a measurement of the latter in the long-baseline neutrino oscillation experiments could 
shed some light on the ball-park magnitude of Mi. 

It is worth mentioning that we have neglected possible renormalization-group running 
effects of neutrino masses and lepton flavor mixing parameters between the scales of v and 
Mi [24]. Such an approximation is expected to be safe in the mi = case, in which three 
light neutrinos have a clear mass hierarchy. In the m 3 = case, in which mi m 2 holds, a 
careful analysis of possible running effects on the FGY ansatz is needed [25] but it is beyond 
the scope of this paper. 

Finally let us comment on possible phenomenological consequences of the FGY ansatz 
for Md in Eq. (4) or equivalently M v in Eq. (6). In the mi = case, one can make use of 
Eqs. (9) and (10) to calculate relevant phase parameters by setting {M v )i2 = 0. A similar 
analysis can be done for the m 3 = case by using Eqs. (14) and (15) and taking (M„) 12 = 0. 
We find that the simple replacements 5 ^ 5 — n and 9 y — > n/2 — 9 y may allow us to write out 
the expressions of a, 0, a, (3 and 7 in the (M„) 12 = case directly from Eqs. (1 1)— (13) and 
(16)— (18). It turns out that the numerical results of a, <f) and a are essentially unchanged, 
but those of (3, 7 and Jcp require the replacements (3 # 7 and Jcp — > — Jcp- The results 
for Y B are essentially identical in (M u ) 12 = and (M v ) 13 = cases. 

IV. SUMMARY 

We have analyzed the minimal seesaw model for leptogenesis and neutrino mixing, in 
which the FGY ansatz is incorporated. We point out a very striking feature of this model: 
all nontrivial CP-violating phases can be determined in terms of the lepton flavor mixing 
angles and the ratio of Amg Un to Am^ tm . This important observation allows us to make 
very specific and testable predictions for the cosmological baryon number asymmetry, the 
effective mass of neutrinoless double beta decay and CP violation in neutrino oscillations. A 
precise measurement of the smallest mixing angle in reactor- and accelerator-based neutrino 
oscillation experiments will be extremely helpful to examine the FGY ansatz and other 
presently viable ansatze of lepton mass matrices. 
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FIGURES 



(a) (b) 




<m> es (eV) sine z 

FIG. 1. Numerical results for the mi = case: (a) dependence of 5, a and on sin^; (b) 
dependence of a, (3 and 7 on sin9 z ; (c) allowed ranges of Jcp and ( Tn^ ee ] (d) dependence of I^b on 
Mi and sin# 2 . The region between two dashed lines in (d) corresponds to the range of 1b allowed 
by current observational data. 
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FIG. 2. Numerical results for the 7713 = case: (a) dependence of 5, a and on sin0 z ; (b) 
dependence of a, (3 and 7 on sin# 2 ; (c) allowed ranges of Jcp and {m) ee ; (d) dependence of Yb on 
M\ and sin^. The region between two dashed lines in (d) corresponds to the range of Yb allowed 
by current observational data. 
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FIG. 3. Numerical illustration of the correlation between Yb and Jcp: (a) in the mi = case 
with Mi = 5 x 10 10 GeV; and (b) in the = case with Mi = 5 x 10 13 GeV. The region between 
two dashed lines in (a) or (b) corresponds to the range of Yq allowed by current observational data. 
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